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Abstract 

We derive the wall crossing formula for the decay of a quarter BPS dyon into a pair of 
half-BPS dyons by analyzing the quantum dynamics of multi-centered black holes in A/" = 4 
supersymmetric string theories. Our analysis encompasses the cases where the final decay 
products are non-primitive dyons. The results are in agreement with the microscopic formula 
for the dyon spectrum in the special case of heterotic string theory on T^. 
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Much of the study of string theoretic black holes has focussed on BPS black holes since the 
degeneracy - or more precisely an appropriate index - associated with BPS states are protected 
and do not change under a continuous variation of the various moduli. Nevertheless it has been 
known for some time [ll|2l|3] that in some cases the index can jump across codimension one 
subspaces of the moduli space on which the original state becomes marginally unstable against 
possible decay into a pair of BPS states. Such subspaces are known as walls of marginal 
stability, and the jump in the index across these walls - known as the wall crossing formula 
- has been the subject of intense investigation in recent years [ll[5l[6l[71[8l[9l[T0l[TTl[T2]. While 
most of the work has been focussed on half-BPS dyonic black holes in A/" = 2 supersymmetric 
string theories, by now we also have a good understanding of this phenomenon in A/" = 4 
and Af = 8 supersymmetric string theories [13l[Tll[T5l[T6l[T7l[18l[T9l[20]. In fact in the latter 
theories the situation is somewhat better since we also have a good understanding of the 
exact spectrum of BPS dyons in these theories. This allows us to verify the general wall 
crossing formula derived from macroscopic considerations involving multi-centered black hole 
solutions [7j against explicit results obtained from microstate counting [T6llT7j . 

For deriving the wall crossing formula for walls on which the decay products carry primitive 
charge vectors, the original techniques of [7] involving analysis of multi-centered black hole 
solutions are equally applicable to theories with A/" = 2, 4 and 8 supersymmetries in four 
dimensions. In the case of A/" = 2 supersymmetric theories these as well as other techniques 
have been developed which allow us to generalize the wall crossing formula to the cases where 
the final decay products carry non-primitive charge vectors |9],|10l[TT]. However so far these 
techniques have not been generalized to A/" > 2 supersymmetric string theories. 

Recently by examining the exact formula for the dyon partition function in the special 
case of heterotic string theory on T^, Refs. [lHl[19] studied the jump in the dyon spectrum of 
this theory across various walls of marginal stability, including the ones on which the decay 
products are non-primitive, and proposed a general wall crossing formula for such walls. The 
purpose of this paper is to give a macroscopic derivation of this formula from the study of 
multi-centered black holes in a general A/" = 4 supersymmetric string theory. As we shall see, 
this is possible and yields results in perfect agreement with the results derived from microscopic 



the leading contribution, this reinforces our belief that black holes capture not only the leading 

^The formula for the dyon partition function was originally guessed in |18pi9j by imposing various consistency 
conditions including known wall crossing formula for decay into primitive dyons. More recently a proof has 
been suggested in pT] , 




Since the jump in the index across these walls is exponentially small compared to 
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contribution to the statistical entropy but also the exponentially suppressed contributions. 

We begin by reviewing the derivation of the wall crossing formula in A/" = 2 supersymmetric 
string theories. Let us consider the wall of marginal stability associated with the decay of a dyon 
of charge (Q, P) into a pair of primitive dyons of charges {Qi, Pi) and {Q2, P2), with Q and P 
denoting electric and magnetic charges in some basis. In this case the jump in the index across 
this wall of marginal stability can be computed from the following simple argument [HIMMTfSH^ . 
A classical analysis shows that a two centered solution of total charge {Q, P), with one center 
having charges {Qi,Pi) and the other center having charges ((52,-^2), exist on one side of the 
wall and does not exist on the other side [HIS]. Thus the jump in the index can be identified 
as the index associated with the two centered solution. One also finds that as we approach the 
wall, the separation between the two centers approaches infinity. Thus in this limit the index 
can be identified as the product of the index associated with each component, and the index 
associated with the supersjTumetric quantum mechanics describing relative motion between 
the two centers. The latter gives a contribution of (— 1)Qi'-P2~Q2-Pi+i ■ P2 — Q2 ■ Pi\- Thus if 
dhiQii Pi) and dh{Q2, P2) denote the index associated with the decay products, the net change 
in the index will be given by 

(_l)Qi-P2-Q2-Pi+l . P2 - Q2 ■ Pll dhiQu Pi) 4(Q2, P2) . (1) 

When the charge vectors {Qi,Pi) and/or {Q21P2) are non-primitive, the above formula is 
known to undergo non-trivial modification [Ql fTOlfTT] . 

Our goal in this paper will be to follow a similar logic for deriving the change in the index 
across a wall of marginal stability where a quarter BPS state in an A/" = 4 supersymmetric 
string theory decays into a pair of half-BPS states. As we shall see, this procedure can give the 
wall crossing formula non only for decay into a pair of primitive dyons but also for decay into 
a pair of non-primitive dyons. We begin by introducing some notations. For a given charge 
vector (Q, P) we define helicity trace -B„ via the relation 



i?„(Q,P) = lrr(Q,p)((-l)^(2/ir) , (2) 

where Tr(^Q^p^ denotes trace over states with charge {Q,P), F denotes the fermion number 
and h denotes the helicity of the state. Since a quarter BPS state breaks 12 out of the 
16 supersymmetries, it has 12 fermion zero modes and we need at least 6 factors of 2h to 
get a non-zero answer for the supertrace. We shall denote by d{Q,P) the index —Bq{Q,P) 
associated with quarter BPS states. A typical quarter BPS supermultiplet contains states 
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of helicities between if — | and if + |, with the state of hehcity H + coming with a 
degeneracy of (^). It is easy to see that the contribution of these states to is given by 
(— 1)^^^^, with the 1/6! in the normahzation factor of cancelhng a factor of 6! coming 
from sum over s. Thus d{Q, P) effectively counts the number of quarter BPS supermultiplets 
carrying charge {Q, P) weighted by (—1)^^, H being the average hehcity of all the states in the 
supermultiplet. On the other hand half-BPS states break 8 out of the 16 supersymmetries, and 
hence we need 4 factors of 2h to get a non-vanishing supertrace. We denote by dh{Qi, Pi) the 
index B4^{Qi,Pi) associated with the half-BPS decay products. dh{Qi,Pi) counts the number 
of half-BPS supermultiplets weighted by (—1)^^, H being the average helicity of all the states 
in the half-BPS supermultiplet. Our goal is to calculate the change in d{Q, P) across a wall of 
marginal stability in terms of dh{Qi,Pi) associated with the half-BPS decay products. 

First consider the case when both decay products are primitive^] In this case on one side 
of the wall of marginal stability we have a two centered classical black hole solution, with one 
center carrying charge (Qiji^i) and the other center carrying charge (Q2,-P2) fl This solution 
ceases to exist on the other side of the wall; hence the jump in d{Q, P) across this wall can 
be identified as the contribution to —Bq{Q, P) from this two centered solution. Now since the 
wall crossing formula is a change in the index, we expect it to be invariant under a continuous 
change in the moduliQ This allows us to work in a region of the moduli space where one of 
the decay products (say with charge (Qi, Pi)) is heavy and the other (with charge {Q2, P2)) is 
light. In this case we can describe the system as the light particle moving in the background 
of the heavy particle, and ignore the backreaction of the light particle on the dynamics of the 
heavy particle. We now note that since the heavy particle breaks 8 out of 16 supersymmetries, 
in order to get a non-vanishing contribution to the supertrace over the states of the heavy 
particle we must insert a factor of (2/i(i))^ into the trace, being the helicity of the heavy 
particle. Since the light particle moves in the background produced by the heavy particle, it 
only feels 8 of the unbroken supersymmetries. Furthermore since the classical two centered 
solution is quarter BPS, the light particle breaks 4 out of these 8 supersymmetries. As a result 



wish to thank F. Denef for discussion on this case. 

^Since both the final dyons are half-BPS states of the full AA = 4 supersymmetry algebra, they are small 
black holes [241125]. but this does not affect our argument. 

^In Af = 2 supersymmetric string theories the decay products are half-BPS and can themselves decay across 
walls of marginal stability. As a result the wall crossing formula also changes as a function of the moduli. In 
contrast in the J\f — A supersymmetric string theories the decay products are half-BPS states whose index does 
not change as we vary the moduli. Thus we expect the wall crossing formula to be unchanged as we vary the 
moduli. 
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it carries 4 fermion zero modes, and we must insert a factor of {2h{2)Y into the supertrace over 
the hght particle degrees of freedom to get a non- vanishing answer. On the other hand since 
in flat space-time the light particle, being half-BPS, breaks 8 out of 16 supersymmetries, it has 
altogether 8 fermion zero modes in flat space-time. 4 of these zero modes must be lifted, i.e. 
take part in the interaction, in the presence of the heavy particle. In fact these must combine 
with the bosonic modes describing the physical coordinates of the light particle to describe a 
supersymmetric quantum mechanics with 4 supersymmetries since the final configuration has 
four unbroken supersymmetries. If d^ei denotes the number of supersymmetric ground states of 
this quantum mechanical system, weighted by (—1)^, then the index —Bq of the two centered 
dyon system will be given by 

dhiQl,Pl)dhiQ2,P2)drel. (3) 

The normalization and sign factors work out as follows. First of all we note that the coefficient 
of (2/i(i))^(2/;,(2))^ term in (2/;,(i) + 2h(2))^ is (^). This factor of (^) combines with the 1/6! in 
the definition of Bq to give l/4!2!. The 4! now cancels the trace of (— l)^(2/;,(i))^ over the 8 
fermion zero modes carried by the heavy state and 2! cancels the trace of (— l)'^(2/i(2))^ over 
the fermion zero modes carried by the light state, leaving behind a minus sign. This minus 
sign compensates for the minus sign in the relation between d{Q, P) and Bq{Q, P). 

In order to calculate drei we can work in a subspace of the moduli space where the original 
dyon with charge {Q,P) and the decay products carrying charges {Qi,Pi) and {Q2,P2) can 
be regarded as half-BPS dyons of an A/" = 2 subalgebra. In this case we can examine the 
situation from the point of view of A/" = 2 supersymmetry. In the absence of the heavy particle 
the light particle would break 4 of the 8 supersymmetries and hence would carry four fermion 
zero modes. However in the presence of the heavy particle these fermion zero modes will be 
lifted since the heavy particle already breaks 4 of the 8 supersymmetries and hence the light 
particle will not break any further supersymmetry. Thus effectively the 4 fermion zero modes 
of the light particle will become interacting and combine with the bosonic coordinates to give 
a supersymmetric quantum mechanics with four supersymmetries. This must be the same 
interacting quantum mechanical system that we got by analyzing the system from the A/" = 4 
viewpoint, - the effect of truncation to the M = 2 subsector being simply the removal of the 
4 non-interacting fermion zero modes on the light state. Now we can use the already existing 
results in A/" = 2 theory to conclude that this supersymmetric quantum mechnics has a set 
of supersymmetric ground states of Witten index (— i)Qi-^2-Q2-Pi+i . — . p^y Using 
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([3]) we now get the total jump in the index across the wall 

Ad{Q, P) = \Q,.p,-Q,.p,\ dh(Q,, P,) d^(Q,, P,) . (4) 

Note that the final formula is symmetric under the exchange of final state charges even though 
we treated them differently in our analysis. 

This finishes our analysis of the wall crossing formula for decay into a pair of primitive 
dyons. Now consider the case where the light dyon carries a non-primitive charge vector, i.e. 
{Q2, P2) is given by some integer N2 times a a primitive charge vector {Q2/N2, P2/N2). In this 
case besides the two centered configuration described earlier, the system also admits multi- 
centered configurations, with the first center carrying charges {Qi, Pi) and the others carrying 
charges {aiQ2/N2,aiP2/N2), with G 'E, ai > 1, — ^2- AH of these configurations 

cease to exist on the other side of the wall and hence could contribute to the wall crossing 
formula. Since we are working in a region of the moduli space where the dyons carrying 
charges of the form {aiQ2/N2,aiP2/N2) are light, they are also weakly interacting. Thus we 
can regard the states of the full system as tensor products of the states of the heavy particle of 
charge {Qi,Pi) and the states of light particles carrying charges {oiiQ2 / N2, aiP2/ N2) moving 
in the background of the heavy particle. However the correponding trace will vanish [12j since 
the heavy particle carries 8 fermion zero modes requiring insertion of 4 powers of helicity into 
the supertrace and each of the light particles carries 4 fermion zero modes, requiring insertion 
of 2 powers of helicity into trace over the Hilbert space of each particle. Since we only have 
a total of 6 powers of helicity in the definition of Bq we cannot saturate all the fermion zero 
modes.. 

There is however an exception to this rule. The above argument assumes that the full 
Hilbert space is a direct product of the Hilbert spaces of the component dyons. However if 
some of the components are identical then we must ( ant i-) symmetrize the wave-function and 
the full Hilbert space is no longer a direct product of the Hilbert spaces of the component 
dyons. Such configurations could give no n- vanishing contribution to the index f] We begin by 
examining the contribution from a multi-centered dyon configurations with one heavy center 
of charge (Qi, -Pi) and L light centers each carrying charges {Q2/ P2I L), L being a factor of 
N2. Furthermore we take all these L light centers in the same internal quantum stated and in 
the same state of the supersymmetric quantum mechanics describing the motion of the light 



^Similar issues arose in the analysis of [2T] . 

^By internal quantum state of a half-BPS dyon we shall refer to supersymmetry singlet part of the state. We 
tensor this state with the states obtained by quantizing the fermion zero modes to get the full supermultiplet. 
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particle in the heavy particle background. We shall argue later that these are the only types 
of multi-centered configurations which contribute to the index. 

Let us now compute the contribution to the index from such a configuration. We denote 
by hint the contribution to the helicity of one of the L light dyons from the internal quantum 
state. This dyon also gets a contribution to the helicity from the supersymmetric quantum 
mechanics describing the motion of the light particle in the heavy particle background, - this 
is given by [7] 

hrei = \Qi ■P2-Q2- Pil/i^L) - (1/2) - integer . (5) 

Finally quantization of the 4 fermion zero modes produces a 4-fold degenerate state, with 1 
state of helicity —1/2, two states of helicity and one state of helicity 1/2. Let us denote 
these four states by \i) with 1 < ^ < 4 and let hi be the contribution to the helicity of the ith 
state from the fermion zero modes. Thus when we take the tensor product of a fixed internal 
state of the light dyon, a fixed state of the supersymmetric quantum mechanics, and the states 
obtained by quantizing the fermion zero mode, the states can be labelled by the index i, with 
total helicity 

hi = hint + Kel + hi . (6) 

We shall refer to these states as single particle states. Our goal is to consider an L particle 
state, with each of the L particles being in the same internal state and same state of the 
supersymmetric quantum mechanics, and calculate the contribution to Bq from these states. 
For this we need to first identify the statistics of the single particle states described above. 
Naively the particle will be bosonic or fermionic depending on whether hi is integer or half 
integer. However in the contribution to hrei given in the part proportional to IQi ■ -P2 — <52 • 
Pi|/(2L) comes from the angular momentum of the electromagnetic field of the dyon system 
and does not directly contribute to the statistics of the light particle. Thus the particle should 
be regarded as bosonic or fermionic depending on whether 2hint + 2hi + 1 is even or odd. The 
L particle states are then labelled by specifying the occupation number rii of the ith state, 
subject to the conditions that rii = L, rii takes values and 1 if 2hint + 2/ij + 1 is odd and 
rij takes all non-negative integer values if 2hint + 2hi + 1 is even. 
Let us now denote hy ql the quantity: 

9L = -^Tr,{{-lf\2hf) , (7) 

where Tr^ denotes trace over the Hilbert space of L particle states introduced above and 
/i in (iTl) stands for the total contribution to the helicity from the L light dyons. Then the 
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contribution to —Bq from these states will be given by the product of ql and dh{Qi, Pi). Using 
the description of the multiparticle states given above, we get 

^ y: i-ir^^-^'^ (2j2n,h)i . (8) 



9l 



Using ([6]) we can write this as 



ni=L 



9l 



9l 



where 



(9) 



(10) 



is most conveniently evaluated by first calculating the partition function 



{rii} 



and then calculating as the coefficient of the e^^ term in 



1 
2 



d(3' 



'12) 



(3=0 



As described earlier, the sum over Ui is restricted to and 1 for 2hint + 2hi + 1 odd and to all 
non-negative integers for 2hint + 2hi + 1 even. This gives 



f{f^,(3) = l[(l-e 



(_l)2'ii„t+2hi 



(13) 



Thus 



ln/(/i (3) = l^_-^Yhint+2hi _ ^/^+2l3(hi„t+hi+hrei)^ 

i 

oo ^ 



k=l 



oo _ 
h 



(14) 



k=l 
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where in the second step we have expanded ln(l — x) in a. Taylor series expansion in x and in 
the last step we have explicitly carried out the sum over i using the fact that there is one state 
with hi = —1/2, two states with hi = and one state with hi = 1/2. This gives, from (fT2l) 

= i^, (15) 

and hence, from ([S]), ([2]), 

_ ^_l^Ql-P2~Q2-Pl+2hint + l ]^ ^iQ^ 

As already mentioned below ([7j), the net contribution to the index —Bq from the specific 
configurations analyzed above will be given by gidhiQi, Pi)- 

The total contribution to —Bq from the configurations where the light dyon is split into L 
identical centers is obtained by summing gi dh{Qi, Pi) over all the internal states of the dyon of 
charge {Q2/L, P2/ L) and all the supersymmetric ground states of the supersymmetric quantum 
mechanics describing the motion of a single particle of charge {Q2/ P2/L) in the background 
of the heavy particle with charge (Qi, Pi). The former gives a factor of dh{Q2/L, P2/L) while 
the latter gives a factor of \Qi ■ P2 — Q2 ■ Pi\/ L [7]. This gives the net contribution to the index 
from these configurations to be 

(_l)QiP2-Q2Pi+i .P,-Q,.P,\ dh(Q,, P,)dh(Q,/L, P2/L) . (17) 

Note that (— l)2'^int has been absorbed into the definition of dh{Q2/ P2/ L). Finally we must 
sum over all possible values of L since all multi-centered configurations with one center having 
charge (Qi, -Pi) and the other L centers having charges {Q2I L, P2/ L) will disappear across the 
wall of marginal stability on which the original dyon decays into dyons of charge (Qi, Pi) and 
(Q2, P2)- This gives the final formula for the jump in the index to be 

(_l)Qi P2-Q2 P:+l . P2 - Q2 ■ Pi I ^^(^1' Pl)dh{Q2/L, P2/L) . (18) 

L\{Q2,P2) 

One could in principle consider more general configurations where the charge vector {Q2, P2) 
splits into different groups with the members within each group being identical but members 
of different groups being not identical. In this case each group will require an insertion of an h"^ 
factor to saturate its fermion zero modes. Since a factor of h^ is already used up by the heavy 
state we see that we can allow at most one group. Thus the configurations we have analyzed 
above are the only ones which can contribute to Bq. 
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This finishes our analysis of the wall crossing formula when one of the decay products is 
non-primitive. What about the case when both decay products are non-primitive? For this we 
let both {Qi, Pi) and {Q2, P2) be non-primitive and continue to work in the corner of the moduli 
space where the dyon of charge (Qi, Pi) is much heavier than the dyon of charge {Q21 -^2)- Now 
the dyon of charge {QiiPi) can also split into multiple centers producing a non-spherically 
symmetric background for the dyons of charge {Q2/ P2/ L). However by arguments similar 
to the ones given above we can conclude that unless all the centers into which ((5i,-Pi) splits 
are in the same quantum state, the contribution to the index from these configurations will 
vanish. If we use position space basis for these centers - which is the natural basis for heavy 
particles - we see that the different centers into which the dyon of charge (Qi, -Pi) splits must 
coincide in space. Thus it continues to produce a spherically symmetric potential for dyons of 
charge {Q2/ P2/ L) and our previous analysis goes through except for a possible change in 
the overall factor associated with the index of the dyon of charge (Qi, Pi). Thus the jump in 
the index must take the form: 

(_l)Q.-P2-Q2-Pi+l .P2-Q2-Pi\ Yl ^'^(^l' Pl)MQ2/L, P2/L) , (19) 

L\{Q2,P2) 

for some function fh{Qi, Pi)- fh{Qi, Pi) gets a contribution of dh{Qi, Pi) from the single dyon 
state but also possible additional contributions from the multi-dyon states into which the dyon 
of charge {Qi, Pi) may split. To determine the form of the function fh we can go to the corner 
of the moduli space where the dyon of charge {Q2, P2) becomes heavy and the dyon of charge 
{Qi,Pi) becomes light and repeat our analysis. This gives the final form of the jump in the 
index for a general decay where both (Qi,Pi) and {Q21P2) are non-primitive: 

Arf(Q,P) = (_i)Qi-P2-Q2-Pi+i \Q,.p^-Q^.p,\ d,{Qi/Li,Pi/Li) Y dhiQ2/L2,P2/L2). 

(20) 

This agrees with the result of [HI [19] derived from the microscopic formula for the dyon 
partition function [T8[ [T9l[2T] for the special case of heterotic string theory compactified on T^. 
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